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Many random variables are almost linear, in the sense that they can be represented 
approximately as the sum of independent components in random phase. Such variables 
(for example, the surface elevation in a random sea) may have a gaussian distribution in the 
first approximation. However in higher approximations the phases of the different com- 
ponents become correlated, due to nonlinear interactions. The purpose of this paper is to 
show theoretically what is the effect of such nonlinearities on the basic gaussian distribution. 

The modified distribution is derived both for a single variable and for two or more 
related variables (such as the components of slope of a random surface). The results are 
applied in the first place to sea waves, and are compared with observations. However the 
analysis is applicable quite generally to any such nonlinear variables. 

Two further problems are solved for weakly nonlinear variables: the mean number of 
zeros per unit time of a stationary random function T(0 and the distribution of the maximum 
values of f(0- These solutions are essentially generalizations of the well-known results of 
Rice for gaussian variables. 

1. Introduction 

Many random variables occurring in physical problems may be considered as the sum of 
a large number of independent components; thus we write 

f=S «& (i) 

where the a t are constants and the £ { are independent random variables symmetically dis- 
tributed about with variance V tj say. Under certain conditions as N-><*> and each F^O, 
the distribution of f tends to a gaussian distribution with variance X) Vi- 
lli particular it is possible to consider some stationary stochastic processes as the limit of 
a sum such as (1), the £* being components corresponding to particular frequencies cr, or wave 
numbers k t ; there may exist a continuous function E(o), the spectral density of f, such that in 
any small interval (a,a-\-da) the sum of the variances V t is given by 



lim 

iV->co 



( 1ZV, ) = E(a)da+0(da) 2 



Such a representation (equivalent to a stochastic integral) has been widely used in the 
theory of noise in electrical circuits [Rice, 1944 and 1945], in the theory of random sea surfaces, 
microseisms, turbulence, and other physical phenomena [Longuet-Higgins, I960]. The repre- 
sentation gives most satisfactory results when the variable f satisfies a linear differential 
equation and can be shown physically to be the resultant of large number of independent 
contributions. 

In some instances, however, the assumption of linearity is only approximately justified, 
and the variable in question may actually satisfy a nonlinear differential equation whose non- 
linear terms are small but not completely negligible. Such is the case with variables f associ- 
ated with sea waves. The distribution p(£) is then nearly gaussian, but not exactly so. The 
question we wish to discuss is: what is the effect of the nonlinearities on p(f), and how can they 
be calculated in terms of the differential equation satisfied by f ? 
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In a recent paper [Longuet-Higgins, 1963], the representation (1) was generalized in the 
following way. Consider the variable 

f =ai&+a*jfefe+aii*&&&+ • • • (2) 

where the a iy a tj , a ijk , etc., are constants and the £< as before are independent random 
variables with variance V t . (The repetition of any subscript will be taken to imply summa- 
tion.) If the a t are given, then by substitution in the differential equation for f it may be 
possible to determine by successive approximation the values of a tj , a iJk , etc. Assuming that 
these are bounded, then as 2V-^°° the contributions to f from successive terms in the series 
can be expected to be of decreasing order of magnitude. 

To calculate p(f) from (2) it is assumed that p(£) is uniquely determined by the sequence 
of its cumulants, which can be found from (2) to any order required. 

This generalization * of the representation (1) leads to a distribution p(f) which is a generali- 
zation of the gaussian distribution ; in fact it is the gaussian distribution multiplied by a sequence 
of Hermite polynomials, of increasing degree but decreasing order of magnitude. The series 
may be only valid asymptotically, and nonuniformly with f . Nevertheless when applied to 
sea waves the second approximation has been shown in [Longuet-Higgins, 1963] to give 
reasonably good agreement with observation. 

In the present paper what we propose to do is to state without detailed proof the general 
results of Longuet-Higgins [1963], and then to apply the results to the solution of two related 
problems: to determine (a) the mean number of zero -crossings of f per unit time and (b) the 
distribution of the values of f at a maximum. 

2. Distribution of a Single Variable 

We may start by writing down the moments of #(?)• Thus taking mean values in (2) 
one has 



The mean values of all odd order terms vanish, while in the terms of even order only those 
remain in which each £ * is paired with a similar f f . Thus 

{=aiiVi+3a i tjjViV j + . . . 

it being assumed that the a's are symmetric in their suffices. Terms involving f f, f % etc. 
are assumed to become negligible in the limit as N— >°° . Similarly since 

t 2 =(a& i +a ij £ i £ j + . . .)(«*£*+ «* *&&+ • • •) 
we have on taking mean values 

t 2 =aia i Vi+(2a ij a J i+a ii a jj )ViVj + 6a i a ijj ViV j J r . . . 

and so on for higher moments. Now some of the terms in the last expression can be factorized, 

e.g., 

a ii a J3 V l Vj= (auVi) (ajjVj) 

while others, e.g., atjajiViVj, cannot. We call terms which cannot be factorized "irreducible". 
For simplicity the following notation is introduced. a tj . . . i is denoted simply by A{, 
where f is the number of suffices i, j, . . . I; and the sum of all irreducible terms in the product 



i A less direct approach is suggested by Wiener [1958]. It should be noted that in the present paper the individual & are not assumed neces- 
sarily to be gaussian . 
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is denoted by (A P A Q . . . A,). Then we have 

P 

F=S[(AA)+(^)U 8 )i 

P, (7 

etc. However, much simpler than the moments of p(f ) are the cumulants k m , defined by 

J_"_ p(f)e««f=esp [|f (i«)+f? (»*)*+ • • •!• 



2! 



We have 



V 



tf 2 =r 2 -(r) 2 =s (4,^,) 



and in general it can be shown that [Longuet-Higgins, 1963] 



K u y , \A p A q . . . A s ) 



(4) 



where p, g, . . . s run through all values from 1 to °° independently. Thus if we retain terms 
up to the third order in the V t we have 



K 1 =(A 2 ) + (A l ) + (A 6 ) 

K 2 = (4?) + (4!) + (At) +2(4,4,) +2(4,4 5 ) +2(4 2 4 4 ) 

K, = (Af) + 3 (A\A 2 ) +3 (4 2 4 4 ) + 6 (4:4,4s) 

^=(4!)+4(4M3)+6(4 2 4^ 



(5) 



We see that K\ and K 2 are both of order V in general, but when n> 2 the lowest nonvanishing 
term in k u is OiV"^ 1 ). If we define the coefficients 



and generally 



\ 4 =KJKI 

K=KJKS' 2 



(6) 
(7) 



we see that X^ is of order V n/2 l in general. 

The density p(£) is now found by inverting the expression (3): 



p(t) ^t XI exp [it (it) +f? (it)2+ • • •] e ~ ittdt: 



Writing 



Kl' 2 ' 
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(8) 



we find [Longuet-Higgins, 1963] that 

where H n (J) denotes the Hermite polynomial of degree n: 

£i n —j l! " 2 "*" 2! 2 2 ' * '* l ; 

This is the required distribution. In the first approximation X 3 , X 4 , . . . are neglected and the 
distribution is gaussian with mean K\ and variance K 2 . In the next approximation a term 
X 3 £T 3 is included, where H 3 is the cubic polynomial (jP— 3/). In the third approximation the 
quartic polynomial H^ and the sextic polynomial H 6 are both involved. Higher approxima- 
tions can be written down at will. 

Equation (9) will be recognized as essentially similar to Edgeworth's form of the 
Gram-Charlier distribution [Edgeworth, 1906]. 

3. Application to a Stochastic Variable 

Suppose now that the variable f is a stationary random function of the space coordinate 
x and the time t, satisfying a nonlinear partial differential equation (or boundary condition) 
in x and t. How is the distribution p(£) related to the spectral density of f ? 

Let the equation satisfied by f be represented symbolically by 

L(t) + QU) + C(t)+ ... -0 (11) 

where L, Q, C, etc., represent operators that are linear, quadratic, cubic, etc., in f. To solve 
(11) for small perturbations we naturally substitute 

r=er (1) +e 2 r (2) +e 3 r (3) + . . ., (12) 

where e is a small parameter. Writing (11) in the form 

£G-)=-Q(r) -<?«■)- • • • (13) 

and equating coefficients of e, e 2 , e 3 , . . . on the two sides of the equation we have successively 

£(f (1) ) = (14) 

m i2) )=-Q,a a) ) as) 

£(f (3> )=-Q(r a) ,f c2) )-c y a (l) ) (i6) 

etc., where Q(f (1 \ f (2) ) denotes an expression that is bilinear in f (1) and f (2) . If there are wave- 
like solutions to (14) then we can write 

f (1> = Sfl. COS (k n -X-a n t+9 n ) (17) 

71 = 1 

where a n and 6 n are amplitude and phase constants and k n and a n are wave numbers and fre- 
quencies respectively. Generally (14) implies a relation between k n and a n . Thus in the case 
of gravity waves on water we have 

rr n == 9\^n\- 

Now (17) can be written 

r (1) =S { a n cos 6 n cos (k n • x— a n t)+a n sin 6 n sin (k n x—a n t)}. 
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<x n t) \ 
<r n t) f 



n=l,N'. 



If now we suppose that a n cos n , a n sin 6 n are normally and independently distributed then 
we have 

22V' 

where 

f n =a„ cos 6 n , a n = cos (k„ • x— <j n t) 

£N+n=a n sin n , a N + n = sin (k n • 

Tims f (1) is expressed in the same form as in (1). 

To find f (2) we now substitute for f (1) in (15). Since the terms on the right are quadratic 
they can in general be expressed as a series of sum or difference wave numbers, i.e., harmonic 
terms in { (k n ±k m )x— (cr w ±o- m )£} . Hence we can in general find expressions for f (2> . On writing 
the solutions again as products of the original harmonics constituents we find 



Thus f (1) + T (2) has the form of (2), as far as the second-order terms. The evaluation of f (3) , 
etc., proceeds similarly. 

The application to gravity waves on water is given in Longuet-Higgins [1963]. It is 
convenient to use the assumption of homogeneity so that we can consider the waves at the 
special point x—0 and time t=0. That is, we may take 



fl, i=l,2,...JV' | 
|^0, i=N' + l, ...2N' ) 



Then for gravity waves it is found that 

1 



(Xij= < 



j= (BTj+Btj-krkj+ih+kj^kikj), i, j=l, 



1 






{BTi-Btr-hk,), i,j=N' + \, ...N 



N' 



where 



B- t 



(V^_V^) 2 -|k 4 -k,| 



B+ _ ( ■&+ ^kjY{k v kj-kjk,) 
(V£+V^) 2 -|k 4 +k,| ' 



The diagonal terms a it are simply given by 

f k i} i=l,...N', 

i=N'H,.. .N. 
So we find, for example, that, to the second order of approximation 






K 1 =a ii V i =j>2k i V i + 5 (-40^=0 

i=l i = N'+l 



(18) 



as was to be expected; while 



S=W«=|} V t = fJE(k)dk 



(19) 
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and 

N' 

= 6 X} onjViVj 

= 6 ffffK(k, k')E(k)E(k')dkdk' (20) 

where K is a function of k and k'. 

In the special case when the waves travel all in the same direction the above expressions 
become very simple, and it turns out that 

K(k,k f )=min (k t k f ). 

Hence introducing the frequency spectrum F(<r) of f we have 

K 2 = CF(<r)(h (21) 

K 3 =6 ffmin (k, k')F{v)F(c')d<jd<r f 

= 12 f °° I C j F(a)da\ F{<r')d<i' (22) 

since k = a 2 /g. 

Returning to the general case, when the waves are not unidirectional, it can be shown 
that the expression on the right of (22), which we denote by /, is related to K z by the inequalities 

0.44 I<K 3 <IM I. 

Combined with (21) this gives us bounds for the skewness coefficient \^=K 3 /Kl /2 , in terms 
of the frequency spectrum F(a) and irrespective of the directional properties of the spectrum. 

These results have been compared with the observations of Kinsman [1960; Longuet- 
Higgins, 1963]. It is found that the inequality on the skewness is satisfied in most cases of 
observed wave spectra. 

Some observed distributions of wave height have been compared by Kinsman with 
a Gram-Cbarlier distribution based on the measured coefficients of skewness and kurtosis. 
As shown in figure 1, these are a better fit to the observations than the simple uncorrected 
gaussian distribution. Kinsman's suggested distribution differs from the form of the Gram- 
Charlier distribution found in section 2 of this paper, since it does not include a term in H 6 . 
Nevertheless the difference resulting from the terms in H 4 and H 6 is so small in his measure- 
ments that the correction is essentially given by the second term X 3 i? 3 . Hence his figure can 
be used effectively as an assessment of the present theory. 

If one attempts to carry the formal calculation of the moments to third and higher 
approximations in the application to water waves, the calculation breaks down. This is 
because of the occurrence of resonant interactions at the third approximation [Phillips, 1960; 
Hasselmann, 1962; and Benney, 1962], which render some of the third-order terms a ij]c in the 
series (2) slowly dependent on the time t. Thus the present method of calculation is consistent 
only as far as the second approximation in water waves. 

In general, however, there seems no reason why the process should not be carried further 
and the higher-order terms in (2) be evaluated. 
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FIGURE AHI 2.05 

FREQUENCY DISTRIBUTION 
OF WATER LEVEL 

RECORD 017 



WIND SPEED AT 1.22 m 3.88 m/s 0.45 m/s 

WIND DIRECTION AT I 52 m 2 77 °T io °T 

ESTIMATED FETCH 17 00 

AIR TEMPERATURE 29.51 

WATER TEMPERATURE 2 7.55 



Gram - Char I i er 

Gaussian 

O Measured Frequency Distribution 

Skewness 0.175 Kurtosis 0.050 

Gaussian X 2 = 38.999 (p<0.00l) 
Gram-Charlier X 2 = 7.466 (p-0.590 




15 20 2.5 3.0 3.5 



Figure 1. From Kinsman 1960. Comparison of an observed distribution of wave height in sea waves with the 
normal distribution (broken line) and with a Gram-Charlier distribution (solid line). 



4. The Joint Distribution of Two Variables 

The method of section 2 can be extended to the determination of the joint distribution 
of any number of variables of the type 



£=<Xi£i+<Xi££j + a ijk ti££ k + . . . i|=ftfe+/3yf*fe+j8ytfit/f*+ • • 



(23) 



where the a's and P's are known constants and the £* are random variables defined as before. 
For simplicity we state the results for two variables f, rj. 

Let us denote by A p and B Q the terms a tj . . . t and B i} . . m where p and g are the 
numbers of suffices i, j . . . I and i, j . . . m; and let {A V B Q . . .) denote the irreducible part 
of the mean product 

(^■i&ifj! • • • Eix) GBfl£< a €'2 • • • W • • •• 

Then it can be shown [Longuet-Higgins, 1963] that the cumulants K mn in the joint distribution 
of f and rj are given by 



J^-mn — 



S (A Pl . . . A Pm B qi . . . B Qn ) 

pi...p m qi...q n 



(24) 
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the summation extending over all positive integral values of the p t and g ; . In particular we 
have K m0 =K m , the cumulant of p(£), as found in section 2, and as far as the terms in V s we 
have 

K n =(A 1 B 1 )+[(A 1 B 3 ) + (A 2 B 2 ) + (A 3 B 1 )]+[(A l B 5 ) + (A 2 B^)+ . . . +(A S B 1 )}, 

K 2 M(AiB 2 )+2(A 1 A 2 B 1 )]M(A 2 1 Bd+2(A l A 2 B 3 )+2(A 1 A s B 2 ) + (AlB 2 )+2(A 1 A i B 1 )+2(A 2 A 3 Bm, 

K 22 =[2(AlB 1 B i ) + (AlBl) +4(A 1 A 2 B 1 B 2 ) + (A 2 2 Bf) +2(A 1 A 3 Bf)], 

K il =[(A\B s +Z(A\A2B 2 )+HA i AlB 1 )]. 

The joint distribution of f , i) is now given by 

piS, V> 
where 



exp[-|(/ 2 -2p// / +/ /2 )/(l-P 2 ] 
2^(^ 20 K 2-^? 1 ) i 



1+c (X30^30 + 3X 2 li?21 + 3\l2-fi r i2+Xo3^03) + - • • (25) 



./= 






/' = 



;--g 01 
P=X„ 



(26) 



" m * (K£J&) 
and H mn {j,f) p) is a two-dimensional analog of the Hermite polynomial, viz 



H mn exp [-|(/ 2 -2p//'+/' 2 )/(l-p 2 )]=(-l)»+" ~^^ exp 



K/ 2 -2p//'+/' 2 )/(i-p 2 )] 

(27) 



In particular 



(28) 



#io=C/-p/')/(l-P 2 ) 

H 20 =(f- P jrKi-pi)-i 

Hu = tf-Pf)(f-PJ)K1-P 2 )+P 

H 02 =(f-pfy/tt-p 2 )-i. 

The first approximation, in which \ mn (m+??>2) are neglected, is the familiar bivariate gaussian 
distribution, as we would expect. In the second approximation cubic polynomials in /, /' 
must be included, which introduce various types of skewness depending on the coefficients 
X30, X21, X12, and X03. In higher approximations we encounter further terms in a bivariate 
Gram-Charlier series. 

In the paper [Longuet-Higgins, 1963] this distribution was applied to the joint distribution 
of the two components of slope of a random surface. In the present paper we shall apply it to 
two different problems: the mean number of zeros per unit time of a nonlinear random function 
f(f) ; and the distribution of the maximum values of f (t). 

5. Mean Number Of Zeros Per Unit Time 

We now consider the problem of determining the average number of zero-crossings (or 
of maxima and minima, which are zeros of the derivative) per unit time in a wave record. 
We assume that the process is stationary to all orders. 

A general formula for the number of crossings of a level f per unit time t is due to Kac 

[1943]. 
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2vo-)= r P (f, r0IM#«, (29) 

«/ —oo 

where p(f, f?) denotes the joint distribution of f and its time-derivative f ,. Now for this 
distribution we have 

2fio=? =0 
and 



K n =tt-h t ~\(n=o 



for any stationary process f, so that p = \ n = and from (25) 

p(t, U)= «~ h " 2+ '\. f"n» {\ 30 H 3 (f)H (f) + . ..} + .. .] (30) 



27t(jBl 2 0^M)2) 

where 






Now since e~* f ' \f\ is an even function of/', the integral of this function multiplied by //„(/') 
will vanish whenever n is odd. Moreover, when n is even and greater than 

C e_v 'V'|fl'.a , )<V / =2j_" <i y'(-i)»^ 1 i (e-i'-W 

=2H„_ 2 (0) 
on integration by parts. So on substitution in (29) we find 

^-^STOl [ 1+ o {^(/ 3 -3/)+3W}+- • •} (3D 

Thus the number of zeros per unit time, considered as a function of f, has a skewness X 30 similar 
to that of p(f), but a mean value 

2 ^12^20 == 2 K12IK02 

different from the mean of f . 

If we take/=0 in (31) then the third-order terms vanish and we have, to fourth order 

N{0)= t(k! K Q2 ^ \} + k {^oH,(0)+6 22 H 2 (0)-\ 04 H (0) }] 

= T(K 20 K 02 y* L 1+ 8 X40 ~4 X22 ~24 X ° 4 j (32) 

6. Distribution of Heights of Maxima 

The distribution p(£ m&x ) of the heights of the maxima of f is given quite generally by 

P(fmax)= f Ptf, 0, fu)|fll|#«I (33) 

where ^>(f , f t , f n ) denotes the joint distribution of f and its first two derivatives with respect 
to t [Rice, 1944-1945; Cartwright andLonguet-Higgins, 1956]. Now if K ijk denotes the (i,j, k)th 
joint cumulant of f, f ,, f tt , it is clear that K m , K i , K 00 i all vanish, as do K nQ and ifon- Hence 
we have (as in sec. 3) 

1057 



+K m t" 2 +2K m tt")+ 1 ^ (K 3O0 t*+3K 21Q tH'+ ...)+.. Pidtdt'dt". 



Writing now 



I 6 /-TVo02j in — ^002,/ >^ 



we have 
2>(f, f«, ?«) = 



+^{X 3 oo(^) 3 +3X 210 (^) 2 (^)+ ...}+.. .]cfe<fe'<fe" 
where p=if 101/(^200^002) 1/2 . Now in 

I I I p -i(fs+f's'H''s")-Hs 2 +s'2+s"2+2pss'')X(is)r(is'y'{is''y"dsds'ds'' 

the terms in s' are separable from the terms in s and s" and hence as in sections 2 and 5 the 
above expression equals 



1 



(1-P 2 ) 



n H r , r ..(j,j"; p)H r ,(f)e-iy"+^- 2 "ff"^" !) ^-^. 



Hence we find 



2?(f, ?*, £«) = 



3 -i[/' 2 +(/ J -2p//"+/" 2 )/(l-p 2 )] 



(2x) 3 <HK 2W K 020 K 002 ) a (1 - p 2 ) * 

x[l+i{X 3 „oH 30 (/,/"; p) #„(/') +3A 21 „H 20 a./"; P)ff,(/)- 



••} + 



•]■ (34) 



Tofind#(f max ) from (33) we set f t =0=j', multiply by |f n | and integrate with respect to $ tt 
from —oo to 0. Writing 

/=(1-P 2 )*Z ] 

/"=-(l-p 2 )^j 

in the resulting integral we have 



2>(£nuu) = 



'(2 w )^(K 200 ^oc 2 ) 



/; 



2/6 



-Kx 2 4-2p^+?/2) 



X 1+tt { (A3oo^3o+3X 2 i fir2i+3Xio2^ r i2+Xoo3^o3)— 3Xi2o^io— 3X 2i^oi}+ • • • \dy (35) 



where 



H 10 =x+py 
H 01 =px+y 

H 20 =(x+p V y-i 

H n =(x J rpy)(px+y) — p 

H 02 =(px+ y y-i 
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and generally 

To evaluate the integrals we use the following results. Straightforward integration by parts 
gives 

/» 00 

yH n , n e-K*>+%"«+rtdy=(H m , a - 2 ) v=0 e-i*> (37) 

whenever «> 2. To deal with the cases n=\, 0, we note that 

So by repeated application of the operator —I ^ p ^— ): 

H 10 -pH 01 =x(1-p 2 ) = G 1 

n 20 -2pii u +p 2 n 02 =x*(i-p 2 y-(i-p*)=G 2 

^3o-3p// 21 + 3p 2 H 12 -p :! /7o3=^(l-P 2 ) 3 -3x(l-p 2 ) 2 =(?3 
and in general 

//,,o-('|'')p^n- I ,,+ • • • = {i- P 2 )>n n {x{\- P >y>)=G„, 
say. Thus in general, 

Jo Jo 

= J[" [(") pfl._,.,-(5) p 2 tf„- 2 . 2 + . . . +(-l)"- , p"fl r ft.+0„] 

y^ e -lU 2 +2 p ry-\-y 2 ) ( / y 
£(") Pffn-Lo-Q P 2 ^-2..+ • • • +(-l)"-yff M -.] ( .0 

+ <W(x;/>) (38) 

where 

F(x; p)= f " e -i^+2A ri /+^)^ = g-i^(i-P») f e -h*dz. (39) 

Jo J PI 

Similarly, when m=0, 

£" Aoe-^WflipfJQ P Hn-i,i-(£) P 2 H n . 2 . 2 + . . .+#„] ri^Wdj 

»p[|( ri Y 1 )p^- 2 ,o-( ri ^ 1 )p 2 ^ n -3,i+ . . .^_ o e-v+Gn-Jix; p)] 
+[-Qp^. 2 ,,+Qp j F„. s , 1 - . . .~] v _ o e-&+G n [e-i*- P xF(x; P )\ 

=[l Q) p 2 H„_ 2t0 -2 (j) p»H,_ M + • • • +0,] «-**" 

+ p(n(? B _ 1 -x(?„)/' 1 (x;p). (40) 
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Using formulas (7.5) to (7.8) we find altogether 

Kr m aO= [go2o( ^5 02 ~/' ol)1/2] Xre-^-p^(x; p) 
{Zir) Jx 2oo-"-oo2 L 



+ix 3 oo{((l-P 2 )V-(3-9p 2 +5p 4 )^e- il - 2 -(p(l-P 2 )V-6(l-p 2 ) 2 a ; 2 +3p(l-p 2 )) J F'^;p)} 
+|x2oi{p(2-p 2 )x e -^ 2 +((l-p 2 ) 2 x 2 -(l-p 2 ))^(a ; ;p)} 

+i X 102 ^^ 2 +i Xoo 3 pa;e^ 2 -| X 120 { (l-p 2 )xe-^+ ( P - P x+p 3 x)F(x; p)}-\ \ m F(x; p)] (41) 



where 



__/ 



.u 



(1-P 2 ) 1 ' 2 [#200(1-P 2 )] 1/2 
K m 



(^200^002) 



The cumulants can be expressed in terms of the spectrum in the following way. If 



ttt= z yiZi+ytj£i%j+ymZiZ&ic+ • • • 



then we may take 



a,= (l,l,...l;0,0,...0) 

0,= (O,O,. . . 0;— «ri,— ir a , . . . -a*) 

t«=W,^ j .. .<^o,o,...o) 



Also if 



(««) = 
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then 




<P») = 



(— o- } a N+ i t N+i-{-<Tiau) ' ' ' (—o'i ( ^n+i,2N-\~^noci,n) 



((TiCKii — (Tia N +i tN+ i) - • • {(T\ai N — <Tn<Xn+\,2n) • • 



(<T N (X Nj i — <Tl<X2N,N+l) ' ' ' (vn<XnN ^N^2N,2n) 




and 



where 



(y tl )- 



/c 0\ 



\0 D/ 



r (—2o 2 1 a n +2o 2 1 a ff+1 , N+ i) 



[— (o"? + ff^)ai,jv+2ffi(7 W arjv + i, 2atP 



and in D the suffices of a fj are each increased or decreased by N. The nonvanishing cumulants 
can now be expressed as follows: 



K m ==(Al)=a i a i V i = \Edo=m 
Kw=(Bl)=PiPiVi= ^o 2 Edo=m 2 
K m2 =(C 2 )=y i y i V t = \o K Edo—m^ 
K m = (A 1 C 1 )=a t y i V i =-jo 2 Edo= 



-m 2 



Also 



and 



K 300 =S(A 2 1 A 2 ) = Qa t a j a tj V t Vj=6 ((a(o, o')EE' 'dodo' 

2f 03t r=3(£tB 2 )=6W„W=0 

Koo3 =^C 2 1 C 2 )=Qy i ymiV i V } =6 jfo 2 o' 2 {-(o 2 +o' 2 ) a (o, o')+2oo' a*{o, <r') } EE'drtfo 

Xwi=(^JCi) =2a i a j y tj V t V i =2 (( {-(o 2 +o' 2 )a(o, o')-2oo'a*(o, o') } EE' dodo' 

K m =(C 2 A 2 ) = 2y i y j a (J V t V i =2jjo 2 o' 2 a(o,o')EE'dodo' 
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K l20 =(B 2 A 2 )=2^a ij V i V j =2 fT<rV 2 a*(<r, o')EE' dcda' 

Zo2i=(5?C 2 ) =2p l P i y {J V ( Vj=2j(<T<r' { - (a 2 +<r' 2 )a*(a, a') +2aa'a(<T, <r') } EE'dada'. 
One can if necessary retain the next term in K m , K m , K m , and K m ', thus 

K 2 oo^(A 2 ) + (AS)=a i a i V i +2a ij a ij V i V j =JEda + 2Jf[a(a,a')fEE'd<rda' 1 

with similar expressions for K 02 o and K m . (Note that K m = —K 020 to all orders). 
In the first approximation (7.9) gives 



m 2 (wiom 4 — ml)* 
(2x) 3/2 m m 4 



?(fj= M nrr [ e -^-^(x; p)i 



where 

x=f/(l — p 2 )?, p= — m 2 /(m mji 

This is the distribution obtained by Rice [1944-1945] and studied by Cartwright and Longuet- 
Higgins [1956]. 

The remaining terms in (7.9) represent the corrections to this distribution, which are 
order V 1/2 . 
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